Anomalous Behavior of Ru for Catalytic Oxidation: 
A Theoretical Study of the Catalytic Reaction CO + ~0 2 — 
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Recent experiments revealed an anomalous dependence of carbon monoxide oxidation at Ru(0001) on oxygen 
pressure and a particularly high reaction rate. Below we report density functional theory calculations of the 
energetics and reaction pathways of the speculated mechanism. We will show that the exceptionally high rate is 
actuated by a weakly but nevertheless well bound (1 x 1) oxygen adsorbate layer. Furthermore it is found that 
reactions via scattering of gas-phase CO at the oxygen covered surface may play an important role. Our analysis 
reveals, however, that reactions via adsorbed CO molecules (the so-called Langmuir-Hinshelwood mechanism) 
dominate. 
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PACS numbers: 68.35.-p, 82.65.My, 82.65.Jv 

The oxidation of carbon monoxide at transition metal 
surfaces is one of the most extensively studied hetero- 
geneous catalytic reactions (see for example [[]],|]|| and 
references therein). Numerous investigations performed 
under ultra high vacuum (UHV) conditions for many 
different metal surfaces have shown that the reaction 
proceeds via the so-called Langmuir-Hinshelwood (L-H) 
mechanism, which means that the reaction takes place 
between chemisorbed reagents. Typical conditions for 
working catalysts, however, are higher pressure and tem- 
perature. Although it has been demonstrated that for 
a number of systems extrapolation of data over a wide 
pressure range is valid pi, such a conclusion cannot be 
generalized. Recent high gas pressure studies (e.g. at 
about 10 torr), for oxidizing conditions (i.e., at CO/O2 
pressure ratios < 1) HHH reported that CO2 production 
over Ru (0001) is anomalous: 

1) The rate of CO2 production was found to be ex- 
ceptionally high - significantly higher than at any 
other transition metal surface. Interestingly, under 
UHV Ru (0001) is by far the poorest catalyst Q. 

2) The measured kinetic data (activation energy and 
pressure dependencies) were found to be markedly 
different to those for other substrates, and in con- 
trast to the other transition metal catalysts, Pt, Pd, 
Ir, and Rh, highest rates occurred for high concen- 
trations of oxygen at the surface. 

3) Almost no chemisorbed CO was detected during or 
after the reaction. 

It was therefore speculated that the Eley-Rideal (E-R) 
mechanism is operational as opposed to the "usual" L- 
H mechanism. In the E-R mechanism, the reaction oc- 
curs between gas-phase and chemisorbed particles. So 
far E-R mechanisms have only been observed experi- 
mentally for somewhat artificial reactions triggered by 
a beam of atomic hydrogen (or deuterium) [[Tj. To gain 
understanding into the drastically different behavior of 
Ru (0001) for the CO oxidation reaction, we carried out 
density functional theory (DFT) calculations, where the 
main approximations are the supercell approach and the 



employed functional for the exchange-correlation interac- 
tion. For the latter we use the generalized gradient ap- 
proximation (GGA) H which is the best justified treat- 
ment to date. Our study represents the first theoreti- 
cal attempt to follow a heterogenous catalytic reaction 
(molecular and dissociative (atomic) adsorption, surface 
reaction, desorption of products) using DFT-GGA and 
an extended surface. In brief, the details of the theo- 
retical approach are summarized as follows: We use ab 
initio, fully separable, norm-conserving DFT-GGA pseu- 
dopotentials M . The GGA is thus treated in a consistent 
way, from the free atom to the solid surface and the re- 
actants. Relativistic effects are taken into account by 
using spin averaged potentials. The surface calculations 
are performed using a (2 x 2) surface unit cell, a four 
atomic layer Ru slab, and a vacuum region correspond- 
ing to thirteen such layers. We use an energy cut-off of 40 
Ry with three special k-points |1C| in the two-dimensional 
Brillouin zone. Convergence tests for O on Ru (0001) in- 
dicated that this basis set provides a sufficiently accurate 
description [ pT| . The adsorbate structures are created 
on one side of the slab [|l2| where we relax the position of 
all the atoms using a damped molecular dynamics p3[ , 
except for the Ru atoms in the bottom two layers, which 
are kept at their bulk-like positions. Details of the cal- 
culations will be published elsewhere JMl . 

It is well known that under UHV conditions, at room 
temperature, dissociative adsorption of O2 results in a 
saturation coverage of Oo ~ 1/2 corresponding to the 
formation of a (2 x 1) structure pa. Recently, from 
DFT-GGA calculations we predicted that an even higher 
coverage should be stable on the surface, namely, a (1 x 1) 
structure with coverage Oo = 1, where the O atoms 
occupy hep- hollow sites fill]. Subsequently, this struc- 
ture was indeed successfully created under UHV condi- 
tions 111] , where it was concluded that formation of the 
Oo = 1 structure from gas-phase O2 is hindered kinet- 
ically |ll],[l6|], but by offering atomic oxygen (or under 
high pressure conditions), this phase can be attained. 
We also carried out calculations involving higher oxygen 
coverages on the surface, as well as geometries involv- 



ing sub-surface oxygen; these structures were found to 
be unstable and metastable, respectively, with respect to 
gas-phase O2. We therefore do not expect them to play 
an important role for the present investigation. We refer 
to Refs. MJTtJ f° r more details. Because the conditions 
under which the particularly high rates of CO2 formation 
occur involve elevated partial gas pressures (and CO/O2 
ratios < 1), there will be a significant attempt frequency 
to overcome activation barriers for dissociative adsorp- 
tion of 02- Thus, it is likely that during reaction the 
oxygen coverage on the surface approaches one mono- 
layer. We therefore initially assume in our investigation 
of the oxidation of CO at Ru (0001) that the (1 x 1) phase 
covers the surface. 

As mentioned above, from the experiments it had been 
speculated that CO may react from the gas-phase with 
adsorbed oxygen (the E-R mechanism). To investigate 
this possibility we first ask whether CO can adsorb on 
the (1 x l)-O/Ru(0001) surface. The sites considered 
were the on-top and fcc-hollow sites, with respect to the 
Ru (0001) substrate, and a bridge site between two ad- 
sorbed O atoms (compare inset of Fig. 1). For each 
site we calculated the energy as a function of distance 
of the molecule from the surface. In these calculations 
the CO-axis is held perpendicular to the surface with the 
C-end of the molecule closest to the surface. At each 
point we fix the position of the C atom and relax the 
positions of all the O atoms and the top two Ru layers. 
The results are shown in the left panel of Fig. 1, where 
we have also considered the path for CO directly above 
an adsorbed O atom. It can be seen that CO experi- 
ences an energy barrier which starts to build up at about 
2.5 A from the surface for all sites, reflecting a repulsive 
interaction with the O-covered surface. Furthermore, it 
is apparent that the surface potential is significantly cor- 
rugated: considering a constant-total-energy surface as a 
function of the lateral position of the CO, we find that it 
exhibits the lowest energy (but always repulsive) over the 
fcc-hollow site. Thus CO tends to avoid the O adatoms 
but will not form a chemical bond with the metal sub- 
strate. The O-covered surface thus prevents reaction via 
the L-H process leaving the possibility for reaction via 
gas-phase CO with chemisorbed O, i.e. the E-R reac- 
tion. In this respect, for the approach of CO directly 
above an adsorbed O atom (full circles in Fig. 1), we 
find that beyond a critical distance, the repulsive inter- 
action turns into an attractive one, and the CO molecule 
and the adsorbed oxygen atom react to form C02- On 
relaxing the position of the C-atom, the CO2 molecule 
then leaves the surface with a significant energy gain of 
«1.95 eV. The associated energetics are shown in the 
right panel of Fig. 1. It is important here to emphasise 
that we have considered all relevant reaction paths for 
CO at the (1 x l)-0/Ru (0001) surface and that although 
on first consideration, Fig. 1 may appear to suggest that 
the favored reaction pathway for CO2 formation is over 
sites away from the adsorbed O atom, in particular the 
fcc-hollow site, this is not the case: Indeed "slow mov- 



ing" CO molecules with low translational energy will be 
"steered" towards the fcc-hollow sites. These molecules 
will however not achieve reaction due to the sizeable en- 
ergy barrier. Fast CO molecules of high translational 
energy, not susceptable to steering effects, which are in- 
cident at sites away from an adsorbed O atom will also 
not react, but will rather be reflected from the surface. 
Thus, to produce CO2 via this mechanism, the results 
indicate that the molecule must "hit", or get very close 
to, an adsorbed O atom. Interestingly, the calculations 
show that there is a physisorption well for CO, as well 
as for CO2, above the surface (barely visible in Fig. 1). 
The wells are very shallow (ss 0.04 eV) and thus they 
will not play a role. It should be noted, however, that 
the calculated depths are likely to be lower bounds be- 
cause the employed exchange-correlation functional does 
not describe the long-range (van der Waals type) interac- 
tions and the physisorption wells are found at distances 
where the true potential energy is likely to be more at- 
tractive than that given by the DFT-GGA calculation. 

A more detailed understanding of the pathway for re- 
action via scattering of CO is obtained by evaluating an 
appropriate cut through the high-dimensional potential 
energy surface (PES); this cut is defined by two variables: 
the vertical position of the C atom and the vertical po- 
sition of the O adatom below the molecule. In order of 
ease of analysis, the CO-axis is initially held perpendic- 
ular to the surface. The resulting PES is presented in 
Fig. 2, where the coordinate system is shown as the in- 
set. For each point we relaxed all the O atoms (except 
that held fixed at Zq>), and the top two Ru layers. The 
repulsive interaction is again evident as CO nears the 
surface. In response to the approaching molecule, the O 
adatom moves in towards the surface: For example, at a 
distance of Zq = 1.9 A , the O atom is displaced inwards 
by Z = 0.2 A . Thus, the impinging CO molecule "hits" 
a "soft wall" . Reaction to CO2 is achieved via an upward 
movement or "hop" of the O adatom by sw 0.4 A towards 
the CO molecule (corresponding to movement parallel to 
the horizontal axis of Fig. 2) and brings the system to 
the transition state of the reaction marked by the aster- 
isk. In view of the similar masses of O and C, it is likely 
that the impinging CO molecule will impart a significant 
amount of energy to the O adatom, thus stimulating its 
vibrations and facilitating its motion (indicated by the 
oscillations in the dot-dashed curve). The newly formed 
CO2 molecule then finds itself in a particularly unfavor- 
able position and is strongly repelled from the surface 
towards the vacuum region with a large energy gain of 
1.95 eV. In the cut through the PES shown in Fig. 2, the 
energy barrier hindering CO2 formation is «1.6 eV. 

The PES of Fig. 2 corresponds to a constrained situ- 
ation of the surface-CO angle. When this constraint is 
dropped, i.e., when the tilt angle of the CO-axis is allowed 
to relax [n8[, we find that the energy barrier is reduced 
to 1.1 eV, and also that the position of the saddle point 
of the PES occurs closer to the surface (by 0.3 A). At the 
transition state (see Fig. 3), the optimum tilt angle with 



respect to the surface normal is found to be 49° which 
corresponds to a bond angle of 131° for the "C02-like" 
complex. Interestingly, this geometry is very similar to 



that associated with the CO^~ ion |19| and to that pro- 
posed for the "activated complex" for the CO oxidation 
reaction over other transition metal surfaces P0| , 

We have thus identified a likely reaction pathway for 
the E-R mechanism. The activation energy barrier for 
this type of reaction appears to be sizable. However, it is 
very similar to those derived from experimental studies 
of CO oxidation reactions at other surfaces which pro- 
ceed via a L-H process, and also for the measurements at 
Ru (0001) the estimated activation energy is comparable, 
namely 0.85 eV |||]||]. On the basis of the present results 
we predict an energy diagram for the E-R mechanism, 
which is shown in Fig. 4. An estimate of the reaction rate 
gives R = 7.5 x 10 6 exp(-l.l/(fc B T)) s _1 which yields at 
T = 500 K, R = 6 x 10~ 5 C0 2 molecules formed per 
surface Ru atom per second |2lJ] which is about 3 x 10~ 6 
smaller than that observed experimentally la] . This indi- 
cates that this mechanism alone cannot explain the par- 
ticularly high CO2 turnover rate. Nevertheless, the rate 
is only about a factor of 10 -3 less than that for the L-H 
process at Pt or Pd H and with molecular beam experi- 
ments this predicted E-R mechanism and associated en- 
ergetics could possibly be measured for the first time for 
the CO oxidation reaction. 

To understand the high reaction rate reported ex- 
perimentally, we turn to another consideration: CO 
molecules might adsorb at sites at which an oxygen atom 
has been removed (e.g. by the above described E-R re- 
action). Indeed, assuming thermal equilibrium of the 
CO + O2 gas and a mixed CO + O adlayer, the law 
of mass action indicates that about 0.03 % of the sites 
of the (lxl) adlayer will be occupied by CO (we as- 
sumed that the O2 and CO partial pressures are equal, 
the temperature is T = 500 K, and the binding energy 
of CO into an O-vacancy of the (1 x 1) adlayer is calcu- 
lated to be 0.85 eV and the adsorption energy of \ O2 K 
(i.e. an O atom) into a vacancy is 1.20 eV. In reality the 
CO concentration will be even higher because catalysis 
is not ruled by thermodynamic equilibrium but by kinet- 
ics, and we find that CO adsorption into an existing O 
vacancy can proceed basically without hindrance while 
O adsorption (from O2) is hindered by an energy barrier 
(see, e.g. Ref. EH]). Therefore the actual percentage of 
surface sites occupied by CO will be somewhat larger. 
For these CO molecules there is a substantial attempt 
frequency to form a CO2 molecule with neighboring O 
adatoms, now by the L-H mechanism which we expect to 
proceed very efficiently due to the relatively weak bind- 
ing energy of both CO (which we calculate to be about 
half that which it has on the clean surface, - and on the 
surface with O-coverages: O < 0.5) and O atoms in the 
high-coverage (1 x 1) adlayer, as well as the close proxim- 
ity of the constituents. We find that the energy gain on 
CO2 formation (of the surface reaction), via this mecha- 
nism is about 0.66 eV J22J; noticeably smaller than that 



of 1.95 eV (see Fig. 4) but still quite significant if com- 
pared to that of « 0.2 eV at Pt (111) and Pd (111) § as 
determined experimentally. 

In summary, we now have the following picture of CO 
oxidation at Ru(0001): with respect to other transition 
metals, ruthenium binds oxygen particularly strongly. 
Therefore, at low oxygen coverages a Ru catalyst disso- 
ciates O2 efficiently, but (in contrast to e.g. Pd) it holds 
the oxygen (and CO) so strongly that reaction to CO2 is 
disfavored. A good catalyst should actuate this dissocia- 
tion but at the same time should not bind the dissociated 
entities too strongly which gives them good capability 
to diffuse and react. Too strongly bound constituents 
would have little reason to react at all. For oxygen in the 
(1 x l)-0 monolayer, the adsorption energy is significantly 
weaker and thus CO2 formation enhanced. Our results 
indicate that this high coverage oxygen phase enables re- 
action via both scattering of CO (the E-R mechanism) 
and by the L-H mechanism where the former may play 
an important role in initiating the reaction. The high 
rate then develops and is maintained by an efficient L- 
H mechanism. Our theoretical results thus explain the 
anomalous dependence of the reaction on oxygen pres- 
sure, as only under sufficiently high oxygen pressure the 
(1 x 1) layer is attained. 

The different mechanisms identified in our study are 
also likely to play important roles for other catalytic re- 
actions. We hope that the detailed predictions and the 
unusual mechanism outlined above, will be tested by ad- 
ditional experiments. 

We wish to thank Martin Fuchs for his help in creating 
the pseudopotentials. 
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FIG. 1. Energy as a function of distance of the C atom, 
Zq, of the CO and CO2 molecules from the surface for the 
various sites tested. The molecular axes are constrained to be 
perpendicular to the surface. The zero of energy refers to the 
situation where CO is far away from the (1 x l)-0/Ru (0001) 
surface {Zc ~6 A). 
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FIG. 2. Cut through the high-dimensional potential energy 
surface (PES) as a function of the positions of the C atom, 
Zc, and the O adatom, Zq (see inset). The molecular axes are 
constrained to be perpendicular to the surface. Positive en- 
ergies are shown as continuous lines, negative ones as dashed 
lines. The contour-line spacing is 0.6 eV. The dot-dashed line 
indicates a possible reaction pathway. 
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FIG. 3. Transition state geometry identified for the reac- 
tion of gas-phase CO with adsorbed oxygen when the con- 
straint on the molecular axis is relaxed. The large, medium, 
and small circles represent Ru, O, and C atoms, respectively. 
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FIG. 4. Calculated energy diagram 
of CO oxidation at Ru(0001). Note that the depths of the 
physisorption wells are exaggerated for clarity. 



